NASA Technical Memorandum 101414 


Noninteractive Macroscopic Reliability Model 
for Ceramic Matrix Composites With 
Orthotropic Material Symmetry 


»C1II*MB*CI3VE t ACEOSCC5IC 
AiLIAEIilH HCDIL iOR CEBAEIC EAIBIX 
CCHPC5I1ES tllH CEIHCIBOPIC 20K 

SKtltJElBX (KA£A) 9 p 


N89-15H37 

Unclas 

0187920 


Stephen F. Duffy 
Cleveland State University 
Cleveland, Ohio 


and 


Jane M. Manderscheid 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

34th International Gas Turbine and Aeroengine Congress and Exposition 
sponsored by the American Society of Mechanical Engineers 
Toronto, Canada, June 4-8, 1989 




E- 451 2 


NONINTERACTIVE MACROSCOPIC RELIABILITY MODEL FOR CERAMIC MATRIX COMPOSITES 
WITH ORTHOTROPIC MATERIAL SYMMETRY 

Stephen F. Duffy* 

Dept, of Civil Engineering and Engineering Mechanics 
Cleveland State University 
Cleveland, Ohio 44115 

and 

Jane M. Manderscheid 
Structural Integrity Branch 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland. Ohio 44135 


ABSTRACT 

A macroscopic noninteractive reliability model for 
ceramic matrix composites is presented. The model is 
multiaxial and applicable to composites that can be 
characterized as orthotropic. Tensorial invariant 
theory is used to create an integrity basis with invari- 
ants that correspond to physical mechanisms related to 
fracture. This integrity basis is then used to con- 
struct a failure function per unit volume (or area) of 
material. It is assumed that the overall strength of 
the composite is governed by weakest link theory. This 
leads to a Weibull type model similar in nature to the 
principle of independent action (PIA) model for iso- 
tropic monolithic ceramics. An experimental program to 
obtain model parameters is briefly discussed. In addi- 
tion, qualitative features of the model are illustrated 
by presenting reliability surfaces for various model 
parameters. 

INTRODUCTION 

The potential advantages of ceramic matrix compos- 
ites include increased fracture toughness as well as 
creep and corrosion resistance at very high service tem- 
peratures. The primary applications under considera- 
tion are advanced turbine engine components, cutting 
tool bits, heat exchangers and aerospace components 
(specifically those of the national aerospace plane). 
Considering that these composites will be produced from 
nonstrategic materials, it is not surprising that con- 
certed research efforts are underway both in the field 
of materials science to advance processing techniques 
and in the field of engineering mechanics to develop 
design methodologies for these material systems. 

The material system of interest in this paper is 
the whisker-toughened ceramic matrix composite. With 
this media the reliability analysis must account for 
material symmetry imposed by whisker orientation. Duffy 
and Arnold (1989) presented a macroscopic model that 
accounted for the transversely isotropic material symme- 
try often encountered in hot-pressed and injection 
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molded whisker-toughened ceramics. A similar approach 
is used herein to develop a noninteractive reliability 
model for a material with orthotropic symmetry. This 
continuum approach excludes any consideration of the 
microstructural events that involve interactions between 
individual whiskers and the matrix. Other authors have 
addressed fracture of ceramic matrix composites on a 
more local scale. Wetherhold (1989) developed a model 
based on probabilistic principles to compute an 
increased energy absorption during fracture due to 
whisker pull-out. Faber and Evans (1983) have addressed 
the process of crack deflection, and Lange (1970) has 
modeled crack pinning. The latter two approaches are 
founded in deterministic fracture mechanics. Knowing 
that these crack mitigation processes strongly interact, 
it is a seemingly intractable task to experimentally 
detect or analytically predict the sequence of mecha- 
nisms leading to failure. A more feasible approach is 
the continuum based criterion whereby reliability is 
computed in terms of macrovariables. 

As pointed out by Leckie (1981), the difference 
between the materials scientist and the engineer is one 
of scale. He notes that the materials scientist is 
interested in mechanisms of failure at the microstruc- 
tural level and the engineer focuses on this issue at 
the component level. We adopt the engineer's viewpoint 
and present a model of practical utility that macroscop- 
i cal 1 y captures the probabilistic failure phenomenon of 
whisker-toughened ceramics. This point of view implies 
that the material element under consideration is small 
enough to be homogeneous in stress and temperature, yet 
large enough to contain a sufficient number of whiskers 
such that the element is a statistically homogeneous 
continuum. This does not imply that the microscopic 
and macroscopic levels of focus are mutually exclusive. 
Indeed a close relationship must exist between the mate- 
rials scientist and engineer so as to develop better 
failure models . 

NONINTERACTIVE RELIABILITY MODEL 

Here we consider a continuum to be a chain com- 
prised of links connected in series. Therefore, the 
overall strength of the continuum is governed by the 
strength of its weakest link. We further assume that 
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the events leading to failure of an individual link are 
not influenced by any other link in the chain. Defining 
f as the failure of an individual link, then 

f * t AV , (1) 

where AV denotes an increment in volume and ^ is a 
failure function per unit volume of material. Taking 
r as the reliability of an individual link, then 

r - 1 ~ t|r AV . (2) 

If the failure of an individual link is considered a 
statistical event, and we assume these events are inde- 
pendent, then the reliability of the continuum, denoted 
as R is 
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As f is a scalar valued function, it must remain 
form invariant under arbitrary proper orthogonal trans- 
formations. Work by Reiner (1945), Rivlin and Smith 
(1969), Spencer (1971) and others demonstrate that 
through the application of the Cayl ey-Hami 1 ton theorem 
and elementary properties of tensors, a finite set of 
invariants (as opposed to the work of Tsai and Wu 
(1971), where an Infinite number of invariants are 
allowed) known as an integrity basis can be developed. 
Form invariance of ^ is ensured if dependence is taken 
on invariants that constitute the integrity basis, or 
any subset thereof. Adapting the above mentioned work 
to y results in an integrity basis composed of 28 ten- 
sor products. Following arguments similar to Spencer 
(1984) several of these tensor products are equal and 
others are trivial identities such that the final integ- 
rity basis for y contains only the invariants 
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I 4 = aj ajaj j , 

(12) 

I 5 = ajajojkakj , 
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Here f(xp is the failure function per unit volume at 
position xj within the continuum. Lower case Roman 
letter subscripts denote tensor indices with an implied 
range from 1 to 3. Greek letter subscripts are associ- 
ated with products or summations with ranges that are 
explicit in each expression. Adopting an argument used 
by Cassenti (1984), the reliability of the continuum is 
given by the following expression 


R = 
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Ig = bibjcr-p (14) 


and 


l 7 = bjbjajkaki . (15) 

A slightly different set of invariants that corre- 
sponds to physical mechanisms related to fracture can be 
constructed from the above integrity basis. This new 
set of invariants includes 


For orthotropic composites the failure function must 
also reflect the appropriate material symmetry. This 
requi res 

V * j,a 1 .b. ) (5) 

where ai and bi are unit vectors that identify 
local material orientations, and represents the 

Cauchy stress tensor. These orientations are depicted 
in Fig. 1. 

The sense of a^ and bi is immaterial, thus 
their influence is taken through the products ajaj 
and bi b j , i.e. 

» *jj( a. j , a • a^ , b * b j ) . (6) 

Note that a^aj and b^ b j are symmetric second order 
tensors that satisfy the identities 

ai * 1 (7) 

and 

bi bi = 1 (8) 
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Considering a uniformly stressed volume, or in the 
context of Weibull analysis a single link, the invariant 

II corresponds to the magnitude of the stress component 

in the direction of a^ , as shown in Fig. 1. I 2 
corresponds to the shear stress on the face normal to 
aj. Similar interpretations can be made for invariants 


Furthermore, the stress and local preferred directions 
may vary from point to point in the continuum. Thus, 

Eq. ( 6 ) implies that the stress field and unit vector 
fields, i.e. a^j(x^), a^x^) and b^x^), must be spec- 
ified to define t* 


I 3 and I 4 and the direction bj . Invariant I 5 is 
the normal stress in the direction defined by the cross 
product of vectors a-j and b^, i.e. 

d i - e ijk a A • (21) 
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where e^ Is the permutation tensor. As proposed, 
these physical mechanisms are independent of Ip and 
I 3 . Taking 
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ensures y is form invariant. 

It is assumed that compressive stresses associated 

with I], I 3 and I 5 do not contribute to failure so 
that 
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In addition, 
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tensile data from fast fracture of simple bend test 
specimens, often referred to as modulus of rupture (MOR) 
bars. The parameters a] and 3i could be obtained 
from MOR bar tests conducted on specimens machined from 
a billet of hot pressed whisker-toughened ceramics. 

These specimens would be oriented along material direc- 
tion ai , shown as orientation 1 in Fig. 2. MOR bar 
specimens machined with orientation 2 (i.e., along b\) 
could be used to determine <13 and [ 33 , and similarly, 
orientation 3 would be used to determine 0:5 and ( 35 . 
The Weibull parameters associated with shear tractions 
across ai and b\ could be obtained from shear tests 
such as Iosipescu tests (Walrath and Adams, 1983). The 
parameters a 2 and 32 could be determined from this 
type of test using specimens with orientation 1 in 
Fig. 2. Similarly, the final two parameters, <14 and 
34 , could be obtained from shear tests conducted on 
specimens with orientation 2 . 

In comparison to the transversely isotropic model 
(Duffy and Arnold, 1989), for which three sets of 
Weibull parameters are necessary, the orthotropic model 
requires five. However it should be noted here that, as 
with monolithic ceramics, it is quite possible that the 
surface and volume of the material will fail due to dis- 
tinctly different flaw populations. These populations 
differ for several reasons: 

( 1 ) the as-fired condition of the surface may be 
different due to the formation of a reaction layer; 

( 2 ) grinding will change the surface and may impart 
subsurface damage; and 

(3) even if the surface flaw distribution is solely 
due to the intersection of internal flaws with the sur- 
face, the presence of the free surface will reduce the 
applied load necessary for fracture. 

The surface would then have different sets of Weibull 
parameters than the volume and Eq. (4) would also have 
to be independently evaluated over the surface area. 


<I 4 ) = I I 4 I <27) 

for all values of I 2 and I 4 . At this point we assume 
that the stress components identified by the invariants 
above act independently (i.e., a noninteractive theory) 
in producing failure. Following reasoning similar to 
Wetherhold (1983), y takes the form 
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Insertion of Eq. (28) into the volume integration given 
by Eq. (4) along with Eqs. (23) to (27) yields a relia- 
bility model for a three-dimensional state of stress in 
an orthotropic ceramic composite. This model is similar 
in nature, yet different in form, to the PIA theory for 
monolithic ceramics. 

In association with each invariant, the a's cor- 
respond to the Weibull shape parameters and the 3's 
correspond to Weibull scale parameters. A variety of 
test methods could be used to determine these model 
parameters. One approach is to obtain the requisite 


IMPLICATIONS OF THE MODEL 

Subsequent to the determination of the Weibull 
parameters, multiaxial experiments should be conducted 
to assess the model's accuracy. One method is tension 
and/or torsion loadings applied to thin-walled tube 
specimens. The thin-walled tube ensures homogeneous, 
biaxial states of stress. The tubular specimen would 
be highly appropriate for planar applications such as 
heat exchangers where one of the material orientations 
(radial) can be ignored. Unfortunately, at the present 
time no data base exists to estimate model parameters, 
although efforts (Shaw and Bubsey, 1987) are underway to 
accomplish this goal. Thus an assessment of the model 
in comparison to experimental data is reserved for a 
later date, and for the examples that follow, model 
parameters are arbitrarily chosen for the purpose of 
i 1 lustration. 

The calculations for the reliability contours shown 
herein are representative of homogeneously stressed con- 
tinuum elements (or links) of unit volume. For dimen- 
sionless R, the Weibull parameter 3 has units of 
stress-(volume) 1/oc . Consider the material symmetry 
where a^ = (1,0,0) and b] * (0,1,0). The uniaxial 
loading 01 ] * 0 (<^2 = a 13 58 °22 = a 23 = a 33 = 0 ) 
this orientation yields the following invariants for y 
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Hence Eq. (4) becomes 



Taking on * 3] in this expression yields R = 0.3679. 
Therefore, any contour R = 0.3679, in a stress space 
containing an, yields (3] as the intercept along the 
ail-axis. This can be seen in Fig. 3(a) where this 
reliability contour is plotted in the on -022 stress 
space, and 3i = 750. Also note that the intercept 
along the cr22~ ax i s • wh i ch represents a uniaxial loading 
in the direction yields the Weibull parameter 

33 » 500. In a similar fashion the pure shear loading 

a] 3 * 0 (o]] = <j] 2 - G 22 * <?23 = a 33 = ^ would yield 

32 as an intercept along the ai 3 -axis, and 

023 * 0(oji - a 12 = CT 13 = a 22 = a 33 = 0 ) would yield 

34 as an intercept along the o 23 -axis. This can be 

seen in Fig. 3(b) where 32 - 400 and 34 = 300. 

Figure 4(a) depicts the intersection of level sur- 
faces of R with the 011-022 stress plane for the 
previously mentioned material orientation. Here 
oq = 15, 03 = 10.5, and as before, 3i * 750 and 

33 = 500. The three surfaces correspond to R= 0.95, 
0.50 and 0.05. Note that a decrease in the a's 
increases the spacing of the contours which indicates a 
higher scatter in fracture strength. This can be seen 
in Fig. 4(b) where a 3 = 5. The spacing between con- 
tours increases in the 022 direction, however the con- 
tours in the u] ] direction remain unchanged because 
there is no variation in oq or 3 ]. If ct 3 had been 
increased, the spacing between contours would diminish, 
and the corners in the first quadrant would sharpen. In 
general, as the a's increase, eventually the reliabil- 
ity contours would not be distinct from each other and 
they would effectively map a deterministic maximum 
stress failure surface. 

We next consider the effects of off-axis (relative 
to loading) material orientation. Figure 5(a) repre- 
sents level surfaces projected onto the 011 -CT 33 stress 
plane for the material orientation defined by 
ai = (1,0,0) and b] = (0,0,1). Here a] - 15, 

3] = 750, a 3 = 10.5 and 33 = 500, and once again three 
contours of reliability are depicted: R = 0.95, 0.5 

and 0.05. Rotation of the material orientation vectors 
such that a j = ( 1 /72 ,0, 1 /72) and b{ = (-1 /72 ,0, 1 /72) , 
for the conditions identical to those associated with 
Fig. 5(a) results in the contours found in Fig. 5(b). 

Due to the new material orientation, compressive compo- 
nents of the stress tensor a ] j contribute to failure. 
This does not contradict the assumption that compressive 
stresses along the material direction ai (identified by 

1 1 ) , and along the material direction b\ (identified by 

I 3 ) do not contribute to a reduction in reliability. In 
this case on and <733 are no longer coincident with 

the stresses associated with I] and I 3 . Expansion of 
the invariants demonstrates that compressive values of 

ail and <733 contribute to both I 2 and I 4 (the shear 
stresses across the material orientation vectors aj 
and bi , respectively) and hence cause a reduction in 
rel iabi 1 i ty. 

Figure 6 (a) represents level surfaces of reliabil- 
ity projected onto 013-011 stress plane. Here a uni- 
axial compressive stress (on) maintains a reliability 
of unity and the contours are symmetric with respect to 
the cn-axis. This results from the assumptions asso- 
ciated with Eqs. (23), (26) and (27). Rotation of the 
material orientation such that ai = (1/72,0,1/72) and 


bj = (-1/72,0,1/72) results in the contours found in 
Fig. 6 (b). Because the contours are closed, any state 
of stress that has nonzero on or <713 components 
would result in a nonzero probability of failure. 

CONCLUDING REMARKS 

In this paper we have developed a reliability model 
for whisker-toughened ceramic composites with ortho- 
tropic material symmetry. The model was constructed in 
a rational manner using an invariant formulation. Such 
an approach indicates the maximum number and form of 
the stress invariants necessary in defining the failure 
function 4 /. A subset of the integrity basis for y 
was constructed with invariants that correspond to the 
macrovariables assumed to be directly related to frac- 
ture. The whiskers in each continuum element were 
assumed to be distributed in such a manner that three 
orthogonal material orientations could be identified. 
These material orientations are not restricted to be 
the same at each point in the continuum and could vary 
along a family of curves within a component. Thus the 
material is locally orthotropic with respect to the 
local material orientation of each element (or link) of 
the component. This offers flexibility when modeling a 
component by a finite element method. 

Relative to conducting experiments to determine 
model parameters, it was pointed out that while tubular 
specimens yield homogeneous states of stress, not all of 
the parameters could be found. However a combination 
of M0R and shear tests would yield all the parameters. 
Furthermore, several figures were presented that 
depicted reliability contours under different loading 
conditions. These figures point to behavior that 
requires verification in multiaxial experiments. Hence, 
a significant amount of experimental data must be gener- 
ated before a final assessment of the theory can be 
made. 
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FIGURE 2. - ORIENTATIONS OF SPECIMENS FOR MOR BAR AND SHEAR 
TESTS FOR FULL THREE-DIMENSIONAL CHARACTERIZATION OF THE 
FIGURE 1, - PHYSICAL INTERPRETATIONS OF INVARIANTS WE I BULL PARAMETERS. 

CORRESPONDING TO FRACTURE MECHANISMS. 
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(B) SHEAR STRESS SPACE %-0 23 WITH P 2 = 400 AND 
P 4 =300. 

F16URE 3. - RELIABILITY CONTOURS OF AN ORTHOTROPIC 
CONTINUUM OF UNIT VOLUME WITH MATERIAL DIRECTIONS 
a; AND bj COINCIDENT WITH THE 1- AND 2-DIRECTIONS 
RESPECTIVELY. FAILURE FUNCTION, t|), EQUAL TO UNITY 


°ir 
(B) 0 3 * 5. 

FIGURE 4. - FAMILIES OF RELIABILITY CONTOURS DEPICTING 
THE EFFECT OF REDUCING THE HE I BULL SHAPE PARAMETER 0, 
MATERIAL DIRECTIONS ARE THE SAME AS IN FIGURE 3 WITH' 
a 1 = 15, p, - 750 AND P 3 = 500. 
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(A) a; AND bj COINCIDENT WITH THE 1- AND 3- DIRECTIONS 
RESPECTIVELY. 


-1500 



(B) MATERIAL DIRECTIONS aj AND bj ROTATED 45° IN THE 
1-3 PLANE. 


FIGURE 5. - LEVEL SURFACES OF RELIABILITY ILLUSTRATING 
THE RESULT OF OFF-AXIS LOADING IN A NORMAL STRESS 
SPACE. WE I BULL PARAMETERS FOR INVARIANTS U AND I, 
ARE THE SAME AS IN FIGURE 3a, ALONG WITH a 2 = 12, 

P 2 88 *100. = 10 AND = 300. 
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a lr MPa 

(B) MATERIAL DIRECTIONS ROTATED 45° IN THE 1-3 PLANE. 
FIGURE 6. - FAMILIES OF RELIABILITY CONTOURS COMPARING 
ON- AND OFF-AXIS ORIENTATIONS OF THE MATERIAL DIREC- 
TIONS aj AND bj IN THE a^-0,3 STRESS SPACE. THE WEI- 

BULL PARAMETERS FOR I^, AND ARE THE SAME AS 
IN FIGURE 5. 
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